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Abstract. Gapped ground states of quantum spin systems have been referred to in the physics 
hterature as being 'in the same phase' if there exists a family of Hamiltonians H{s), with finite range 
interactions depending continuously on s £ [0, 1], such that for each s, H{s) has a non-vanishing 
gap above its ground state and with the two initial states being the ground states oi H{0) and -ff (1), 
respectively. In this work, we give precise conditions under which any two gapped ground states of 
a given quantum spin system that 'belong to the same phase' are automorphically equivalent and 
show that this equivalence can be implemented as a flow generated by an s-dependent interaction 
which decays faster than any power law (in fact, almost exponentially). The flow is constructed 
using Hastings' 'quasi-adiabatic evolution' technique, of which we give a proof extended to infinite- 
dimensional Hilbert spaces. In addition, we derive a general result about the locality properties of 
the effect of perturbations of the dynamics for quantum systems with a quasi-local structure and 
prove that the flow, which we call the spectral flow, connecting the gapped ground states in the same 
phase, satisfies a Lieb-Robinson bound. As a result, we obtain that, in the thermodynamic limit, 
the spectral flow converges to a co-cycle of automorphisms of the algebra of quasi-local observables 
of the infinite spin system. This proves that the ground state phase structure is preserved along 
the curve of models H{s), < s < 1. 



1. Introduction 

Since the discovery of the fractional quantum Hah effect [50] and its description in terms of 
model wave functions with special 'topological' properties [31], there has been great interest in 
quantum phase transition [38] . Experimental and theoretical discoveries of exotic states in strongly 
correlated systems |13j and, more recently, the possibility of using topologically ordered quantum 
phases for quantum information computation [30], have further increased our need to understand 
the nature of quantum phase transitions, and especially of gapped ground states. It is natural to 
ask whether gapped quantum phases and the transitions between them can be classified. The first 
and simplest question is to define precisely what it means for two gapped ground states to belong to 
the same phase. A pragmatic definition that has recently been considered in the literature declares 
two gapped ground states of a quantum spin system to belong to the same phase if there exists a 
family of Hamiltonians H{s), with finite range interactions depending continuously on s € [0,1], 
such that for each s, H{s) has a non-vanishing gap above its ground state, and the two given states 
are the ground states of H{0) and H{1). In other words there is a family of Hamiltonians with 
gapped ground states that interpolate between the given two [lOl [11]. In this paper we prove a 
result that supports this definition. We show that any two gapped ground states in the same phase 
according to this definition are unitarily equivalent, with a unitary that can be obtained as the flow 
of an s-dependent quasi- local interaction which decays almost exponentially fast. When applied 
to models on a finite-dimentsional lattice, this quasi-local structure is sufficient to prove that the 
unitary equivalence of finite volume leads to automorphic equivalence at the level of the C*-algebra 
of quasi- local observables in the thermodynamic limit. 

In statistical mechanics, lattice models with short-range interactions play a central role. Many 
examples of Hamiltonians that can be considered as a perturbation of a model with a known ground 
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state that is sufficiently simple (typically given by finite number of classical spin configurations), 
have been studied by series expansion methods [271 E 121 EH [Ml [291 [SI [IH [13 [51]. Perturbation 
expansions provide a detailed understanding of the ground state under quite general conditions. Of- 
ten, one can prove the existence of a finite correlation length and a non-vanishing spectral gap above 
the ground state, and short-range correlations can in principle be calculated to arbitrary precision. 
The perturbation series one employs in such situations has the structure of a cluster expansion in 
which the cluster geometry is based on the underlying lattice structure and the short-range nature 
of the interactions. The effect of the perturbations can then be understood as approximately local 
modifications of the ground state of the unperturbed model. 

In this work, we start from a different perspective. Suppose we have a family of models defined 
in terms of an interaction <^(s) which depends on a parameter s G [0, 1]: 



Here, A is a finite subset of the lattice T (e.g., T = Z*^), and $(s) is a short-range interaction 
depending smoothly on a parameter s (see Section [5] for the precise conditions on the decay of the 
interactions that we assume). 

Suppose that for all s the ground state(s) of this family of models are isolated from the rest 
of the spectrum by a gap. We prove that the ground state(s) of H\{s) can be obtained from the 
ground state(s) of H\{0) by a unitary transformation Ua{s) which has a quasi-local structure in 
the sense that U\{s) can be regarded as the flow generated by a quasi-local parameter-dependent 
interaction ^{s) which we construct. In the works cited above the goal was to develop a suitable 
perturbation theory which would allow one to prove the existence of a non-vanishing spectral gap, 
among other things. To do that one has to start from a sufficiently simple model at s = and also 
assume that the perturbation potential is sufficiently small. Note that no such smallness condition 
is required on $(s) here. We now make some comments on the methods used in this paper. 

In his 2004 paper |19j Hastings introduced a new technique, which he called 'quasi-adiabatic 
continuation' (see also [M])- He used it in combination with the propagation bounds for quantum 
lattice dynamics originally due to Lieb and Robinson [32] to construct and analyze the variational 
states needed for the proof of a multi-dimensional version of the Lieb-Schultz-Mattis theorem |33j . 
The quasi-adiabatic continuation technique was subsequently elaborated upon and used in new 
applications by Hastings and collaborators [HI [201 [231 [HJ [7] ^ ^^11 as other authors [lOllll]. In this 
paper we give a general account of this technique and show how it allows one to exploit locality 
properties of the dynamics of extended quantum systems with short-range interactions without 
resorting to cluster expansions. 

The starting point of the analysis in all the works mentioned above is a version of our Proposition 
12. 4[ This result shows that the spectral projection associated with an isolated part of the spectrum 
of a family of self-adjoint operators H{s) depending smoothly on a parameter s, can be obtained 
as a unitary evolution. Concretely, let /(s) C M denote an interval such that for all s the spectrum 
of H{s) contained in I{s) is separated by a uniform gap 7 > from the rest of the spectrum of 
H{s), then there exists a curve of unitary operators U{s) such that 



When we apply this result to families of Hamiltonians H\{s) of the form (II. Ih . i.e., with a quasi- 
local structure, we find that the unitaries U\{s) then have the structure of a quasi-local dynamics 
itself. Explicitly, 



(1.1) 




XcA 



P{s) = U{s)P{0)U{s)* 



— C/a(s) = iDA{s)UA{s) 
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where -Da(s) is a self-adjoint operator with the structure of a time-dependent Hamiltonian,z.e. , 
there is an interaction ^{s) such that 

xga 

Because of this quasi-local structure, the flow on the algebra of observables defined by conjugation 
with the unitaries U/^{s), i.e. : 

(1.2) af{A) = UX{s)AUa{s) 

satisfies a propagation bound of Lieb- Robinson type (see Section U]). These propagation bounds — as 
a second application of Lieb-Robinson bounds in this paper — can be used to prove the existence 
of the thermodynamic limit (Section [5]). The main result of this paper is Theorem 15.51 Stated in 
words, it says that if for a differentiable curve of Hamiltonians of the form (jl.ip the gap above the 
ground states does not close along the curve, then, for each s there is an automorphism as of the 
algebra of quasi-local observables which maps the ground states at s = to the ground states at 
s. In particular the simplex of infinite- volume ground states for all values of s is isomorphic to the 
one for s = 0. 

We find the designation 'quasi-adiabatic' of the flow somewhat misleading since there is 
nothing adiabatic about it. The flow does, however, follow the spectral subspace belonging to the 
isolated interval I{s). We will therefore call it the spectral flow. 

This paper is organized as follows. In Section[2l we give a rigorous and self-contained presentation 
of the construction of the spectral fiow in a form that allows for applications with an infinite- 
dimensional Hilbert space. A number of applications where the infinite-dimensional context has 
proven useful have already been considered in the literature, see e.g. [121 ES [Ml [31 [Ml [46]. We 
expect that more applications will be found. 

In Section [3] we use Lieb-Robinson bounds to obtain a locality property of the spectral fiow and 
prove that local perturbations perturb locally in the sense that the dependence of gapped ground 
states (or any other isolated eigenstates) on any given local term in the Hamiltonian is significant 
only in a neighborhood of the support of that term. In Lemma 13.21 we generalize the notion of 
normalized partial trace to infinite-dimensional Hilbert spaces. 

In the final two sections we consider quantum lattice models, or more generally, models defined on 
a metric graph (satisfying suitable conditions) with sufficiently fast decaying interactions. Section 
[H is devoted to showing that the spectral fiow can be generated by time-dependent Hamiltonians 
defined in terms of local interactions. As a consequence, this fiow then also satisfies a Lieb-Robinson 
bound. In Section \5\ we restrict our attention to quantum spin systems, and use the results of 
Section [H to obtain the existence of the thermodynamic limit of the spectral fiow as automorphisms 
on the algebra of quasi-local observables. We conclude the paper with a brief discussion of the 
notion of 'gapped ground state phase', which has been a topic of particular interest in the recent 
literature. 

2. The curve of spectral projections for an isolated part of the spectrum of a 

Hamiltonian with a parameter 

We consider a smooth family of self-adjoint Hamiltonians H{s) = H{s)* parametrized by 
s € [0,1], acting on a Hilbert space Ti. We do not assume that H{s) itself is bounded but the 
s— dependent portion should be. We are interested in the spectral projection P{s) associated with 
an isolated part of the spectrum of H{s). Explicitly, our main assumption on H{s) is the following. 

Assumption 2.1. H{s) is a densely defined self-adjoint operator with bounded derivative H'{s), 
such that ||i/'(s)|| is uniformly bounded for s € [0, 1]. Furthermore, we assume that the spectrum, 
S(s) of H{s) can be decomposed in two parts: S(s) = Si(s) U S2(s); such that inf{|Ai — A2I | Ai € 
El, A2 G S2} = 7 for a constant 7 > 0, uniformly in s. We also assume there are compact intervals 
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I{s), with end points depending smoothly on s and such that Ili(s) C /(s) C (M \ S2(s)), in such 
a way that the distance between I{s) and ^2(3) has a strictly positive lower bound uniformly in s. 

Typically, we have in mind a family of Hamiltonians of the form H{s) = H{0) + ^{s), with 
H'{s) = ^'{s) bounded. Specifically, if H{s) is unbounded, this is due to H{0), which is obviously 
independent of s. Let Ex{s) be the spectral family associated with H{s) and let P{s) := Jj^^^ dEx{s) 
be the spectral projection on the isolated part of the spectrum Si(s). 

The formulation of the main result of this section uses a function w-y G L^(M), depending on a 
parameter 7 > 0, with the following properties. 

Assumption 2.2. G L-^(M) satisfies 

i. is real- valued and / dt'w.y{t) = 1, 

ii. The Fourier transform is supported in the interval [—7,7], i.e., w-y{uj) = 0, if |a;| > 7. 

Such functions exist and were already considered in [21]. In the following lemma, we present a 
family of such functions derived from [251 [IB] and give explicit bounds on their decay that we will 
need in this work and which may also prove useful in future applications. 

Lemma 2.3. Let 7 > and define a positive sequence (an)n>i by setting an = ai(nln^n)~^ for 
n > 2, and choosing ai so that Yl'?^=i = 7/2. Then, the infinite product 

(2.1) =.,,^(!i;i^)^ 

\ tint / 

n=l 

defines an even, non-negative function € L^(M), and we can choose such that J w-y{t)dt = 1. 
With this choice, the following estimate holds. For all t > e^/^/7, 

(2.2) o<^«^(t) <2(e7)'t-exp(^-^^^^ 

Proof. Without loss of generality, we shall assume t > 0. Since each term of the product lies 
between and 1, and by Stirling's formula, 
N . 

w^it) <c,ll ^ c-,(iV!)2ln4^(Ar)(ait)-2^ < Inc.NN^'' In^^ {N)iaity^^e~''' . 

n=l 

The desired bound is obtained by choosing N = [oit/ ln^(7t)J and noting that 7/7 < a\ < 7/2 
and 7/(27r) < < 7/7r. The bounds on a\ follow directly, while the latter estimates are proven 
e.g. in 0]. For t > ei/^/7, TV < 7* so that 

2 7t 
"7ln2(7t), 

Finally, this decay estimate and the a priori bound w^{t) < 1 for all t imply that G L^(M). 

□ 



w^{t) < 2(e7) t ■ exp 



Since the Fourier transform of sin(aa;)/(aa;) is the indicator function of the interval [—a, a], the 
support of corresponds to [—25", 25], where S = fln, and thus (ii) of Assumption 2.2 also 

holds. Moreover, this lemma shows that the function can be chosen to decay faster than any 
power as t —7- 00. This will be important for some of our applications. We can now state and prove 
the main result of this section. 

Proposition 2.4. Let H{s) he a family of self-adjoint operators satisfying Assumption \2. 1\ Then, 
there is a norm- continuous family of unitaries U{s) such that the spectral projections P{s) associated 
with the isolated portion of the spectrum Si(s), are given by 

(2.3) P{s) = U{s)P{0)U{sy. 
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The unitaries are the unique solution of the linear differential equation 
(2.4) -iA.u{s) = Dis)U{s), UiO) = l, 

where 

/oo rt 
dtw^it) / due^"-^(^)i/'(s)e-'"^(^). 
oo Jo 

for any function satisfying Assumption \2.SX 



It is obvious from (j2.5p and the assumption that Wry(t) G M, that D(s) is bounded, self adjoint, 
and the equations (j2.3|) and (j2.4p can be combined into 



(2.6) ^P(,) = i[Z)(.),P(,)]. 

Moreover, boundedness of D(s) implies that the unitaries U (s) are norm continuous. 

The existence of a (bounded holomorphic) transformation function V{s) such that P{s) = 
V{s)P{0)V{s)~^ is a direct consequence of the smoothness of -P(s), see e.g. [26]. The interest 
of the proposition stems from having an explicit formula of a unitary family U{s), from which 
interesting properties can be derived. This constructive aspect is essential for the applications we 
have in mind (see Sections [3] and U]). 

Proof. On the one hand, 

(2.7) P(s) = --L/' dzR{z,s), 

Jr{s) 

where R{z, s) = {H[s) — z)~^ is the resolvent of H{s) at z, and the contour r(s) encircles the real 
interval I{s) in the complex plane. Therefore, 

(2.8) P'{s) = — —[ dzR'{z,s) = —[ dzR{z,s)H'{s)R{z,s), 

27ri Jr(s) 27ri yr(s) 

where the first equality follows by noting that the smooth dependence of s i— )• I{s) and the uniform 
lower bound on the gap imply that the contour T{s) can be kept fixed while differentiating; Namely 
for e small enough, T{s) can be chosen so that it encircles all intervals /(c), a € [s,s + e]. The s 
dependence of T can therefore be taken as purely parametric. Since P{s) is an orthogonal projection, 
P{s)P'{s)P{s) = (1 - P{s))P'{s){l - P{s)) = and therefore. 



P'(s) = — f dz (P(s)R(z, s)H'(s)R(z, s)(l - P(s)) + (1 - P(s))R(z, s)H'(s)R(z, s)P(s)) 
2vn Jr{s) 

(2.9) =r^/ dz [ di^l dX^—-^{dE^{s)H'{s)dEx{s) + dEx{s)H'{s)dE^{s)) 

2vri7r{s) Ji{s) Jm/i(s) fi-zX-z^ 

(2.10) =-/ d^l^ dX-^{dE^is)H'is)dExis)+dEx{s)H'{s)dE^{s)) . 

Ji{s) Jr/i{s) a - /i 

In order to justify the last equality, we interpret the double spectral integral as a double operator 
integral, see e.g. [5], Theorem A.l(iii). Eq. (12. 9p corresponds to the factorization of the symbol 
(/>(A,//) = {X-^l)-'^ of (l2lOD,the auxiliary measure space being {S^,d-i{t)) where S'^ 3 t ^ 7(t) G C 
is a parametrization of r(s). The uniform integrability conditions are met because of the finite size 
of the gap. On the other hand, 

ms).P{s)] = i((l - P{s))D{s)P{s) - P{s)D[s){l - P[s))) 
(2.11) 

= i / diij dx[dtw^{t)[ du{e'''^^-''UExis)H'{s)dEf,{s) -e-'''^^-''UE^{s)H'{s)dExis)) 

Jl(s) Jr/I(s) J Jo 
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which yields (j2.10p after the time integrations are performed, namely 

i J dtw^iit) j\ue^'''^^-i''^ = ±J dtw^{t)j^^{e^'*^^-^^ - 1) 



(2.12) W^{t) :-- 



A — /i ' A — 

where we used first that J tM^(t) = 1 and then the compact support property of UU/ together with 
the fact that |A — /i| > 7 by Assumption (j2.ip . □ 

We now introduce the weight function 

^ ft"^ dCw^fiC) t>0 

which will play a central role in the following applications. As w-y G L^(M), Wj is well-defined. 
Lemma 2.5. For a > define 

-a-O- 
UaiV) = e 1" " , 

on the domain r] > 1. For all integers k > and for all t > such that also 

t 



we have the bound 



>2k + 2, 



^v'uMdv<^^^t''^'uait). 
a 



Proof. For r/ > e^, the function 

in r/ 

is positive, differentiable, and monotone increasing, and 

dr] 11 It? t] \ rj 



dr all 



< - 



ln(77) / 

If we further require 77 > e^, we can also use the bound 1 < 77/(logr7)^, and therefore 



> 2 2 
r] \ T 



By making the substitution to the integration variable r in the integral, we find 



/oo ^ 
v'Mv)dv < ^2^r(2A: + 3,T(t)) 



where the incomplete Gamma function T{n + l,x) can be computed for any integer n > by 
repeated integration by parts: 



/•OO K 

T{n + l,x)= / T"e-^dr = n!e-^ J^l^. 

■^^ k=o ' 

For X > n, this yields the bound 

r(n + l,x) < (n + l)x"e"'^, 
which can be applied with n = 2A; + 2 and x = r(i) < at to conclude the proof. □ 

Lemma 2.6. Let 7 > and the function defined in 12. Then eq. i2.12\) defines a bounded, 
odd function G L-'^(M) with the following properties: 
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i. is continuous and monotone decreasing for t >0. In particular 
(2.13) \\W4oo = W^{0) = l/2; 

ii. \W^{t)\ < G(^)(7|t|), with G(^)(?7) defined for t]>0 by 



(2.14) G(^)(r/) 



i 0<ri<ri* 
35e^r/^n2/7(?7) rj > rj* 



where rj* is the largest real solution of 

35e%^n2/7(??) = 1/2 . 

iii. There is a constant K such that 

(2.15) \\WJ\i<—. 

7 

iv. For t > 0, let 

/CO 
d^W^iC). 

Then, \I^{t)\ < G(^)(7|t|), where G^^HC) is defined for ( > by 



1 I f < c < C* 



2 



7 [l30e2CiV/7(C) OC 
with K as in (iii) and a C* > 0. 

Remark 2.7. It is straightforward to estimate the values of the constants Q* , and by numerical 
integration. One finds 14250 < ?]* < 14251, 36057 < C* < 36058, and K ~ 14708. 

Proof, i. > 0, even, and J w-^/ = 1. With the definition (j2.12p of W^, this imphes 

roo poo I 

(2.17) \W^it)\ < / w^iO d^ < / ^^^(0 = W^{0) = -. 

J\t\ Jo 2 

ii. The bound ()2.2p for gives 

foo roo foo 

\W^{t)\ = / d^w^iO < 26^2 / dau2/7{lO = 2e2 / d?7??^^2/7(^) • 

J\t\ J\t\ J^\t\ 

With k = 1 and a = 2/7, the conditions of Lemma 12.51 are satisfied for 'y\t\ > 561, so that 

(2.18) \W^it)\ < 35e\^\t\)\2Mj\t\) , if j\t\ > 561. 

Using the decay of Uair]) for ry > e^ and the fact that the RHS of ()2.18p exceeds the a priori 
bound (j2.13p for 7|t| = 561, the result follows. 

iii. By (ii) Wi € L^(]R) and |VF7(t)| < |Wi(7t)|, which imphes the existence of a constant K as 
claimed. Using the oddness of and the exphcit function G^^\i]), we choose 

/•oo 

K = ri* + 70e2 / rj'^u^pir]) dr]. 

Jt]* 

iv. Follows by (iii) and another application of Lemma 12. 5[ 

□ 

A straightforward corollary of the decay conditions of the weight function is the following equiv- 
alent form of the generator D(s), eq. (12. 5p . 
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Corollary 2.8. The conclusions of Proposition \2.4\l hold with 

/oo 
-oo 



;Wy{t) = -W^{t) , 



with W-y as in lemma [2. 

Proof. This follows by a simple integration by parts from (|2.5p . By definition of the function W-y , 
we have, for any t G R \ {0}, 

di 

which can be extended by continuity at t = 0. Proposition (j2.4p then yields 

D{s) = - W^{t) / dne^"^(^)F'(s)e-^"^(^) + / dtW^{t) ■ (^'"^'^ H' {s)e-''"^'^ . 

JO — oo J — oo 

The boundary term vanishes by Assumption (j2.ip and the decay of W^. □ 

3. Local perturbations 

The aim of this section is to combine the evolution formula of Section [2] with Lieb- Robinson 
bounds to show that the effect of perturbations with a finite support X can be, to arbitrarily 
good approximation, expressed by the action of a local operator with a support that is a moderate 
enlargement of X. In principle, the following lemma suffices to turn Lieb-Robinson bounds into an 
estimate for the support of a time-evolved observable. 

Lemma 3.1 ([13]). Let Hi and ^2 Hilbert spaces and suppose e > and A € B{'Hi (8) ^^2) are 
such that 

\\[A,t(^B]\\ < e\\B\\ for all B e B{'H2). 
Then, there exists A' € B{'Hi), such that 

(3.1) \\A'®t-A\\<e. 
If dim?^2 < 00, one can simply take 

A' = ^ Tr-w,^, 

dim 7^2 

as is done in [9l[36] (or see (i) in the proof of Lemma 13.21 below) . 

For the applications we have in mind, we want the map A^ A' to be continuous in the weak op- 
erator topology. In finite dimensions the partial trace is of course continuous. In infinite dimensions 
we cannot use the partial trace and the continuity is not obvious. Moreover, it will be convenient 
for us to have a map A' = 11 (A) that is compatible with the tensor product structure of the algebra 
of local observables of a lattice system (see Section [4.ip . For this purpose, we fix a normal state 
p on B{U2) and define the map H : B{ni) ® B{H2) ^ B{Hi) ^ B{ni) ® 1 C B{ni) ® ^(^^2) by 
n = id (81 p. Although the map 11 depends on p, we have the following estimate independent of p. 

Lemma 3.2. Let Hi and T-L2 be Hilbert spaces and suppose e > and A G B{7ii ® ^^2) o,re such 
that 

\\[A,10B]\\ < e\\B\\ for all B e B{n2). 

Then, 

(3.2) ||n(A) - A\\ < 2e. 
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Proof, (i) First, assume dim 7^2 < oo. Then it suffices to take for A' the normalized partial trace 
of A: 

Note that 

A'(^l= [ dU il0U*)A{l®U) 

Ju(H2) 

where dU is the Haar measure on the unitary group, U{'H2)^ of 7^2- Then, by the assumptions of 
the Lemma, one has 

\\A'®t-A\\< [ dU\\{l'S)U*)[A,{l(E)U)]\\ <e. 

(ii) In the case of infinite-dimensional 1-12, we start by defining, for 7] G 7i2, \\ri\\ = 1, € B{T-Li) 
by the formula 

{(t),Ar^ij) = {(l)®r],A'>p®ri), (/>,Ve?^i. 
For 7?,^ G 7^2, let \Cj{ri\ denote the rank-1 operator defined by |^)(?7| 4> = for all </> € 7^2- For 

any three ri,^,x G T~(-2, = ||^|| = ||x|| = 1) note that 

(3.3) A^ ® \r,){x\ = (10 |r?)(e|)^(]l ® \0{x\) ■ 

This equation is easily verified by equating matrix elements with arbitrary tensor product vectors 
(j)® a and ij) ® 13. By the assumptions we then have 

11(1 m\) [A, 1 \iM] (1 ® m\)\\ < II [A, 1 io(??i] II < e . 

By expanding the commutator and simplifying the products in the left hand side of this inequality 
and using (I3.3p we obtain 

(3.4) \\A^-A4 = \\A^®\7^m-Ar,(^\l^m\\<e. 

Next, consider finite-dimensional orthogonal projections P on 7^2- Since, for each such P, 

P)A{1(^ P),l® {PBP)]\\ = \\[{1®P)[A,PBP]{1®P)\\ < \\[A,PBP\\\ < e\\Bl 
by (i), there exists Ap € BiHi) such that 

(3.5) \\Ap(S)P-{liS)P)A{l'S)P)\\<e. 

Explicitly, if xi^ • • • jXn is an o.n. basis of ran P, the construction in part (i) provides 

1 " 

^P = -Y.^x,^ andPp||<p||. 

k=l 

The diameter of the convex hull of {A^^ \ x £ "^2, llxll = 1} is bounded by e due to (|3.4p . It follows 
that for any two finite-dimensonial projections P, Q on 7i2 

ll^p -^qII < e- 

Now, we prove the bound: 

\\Ap ®1-A\\<2e 

by contradiction. Suppose that for some P, \\Ap ®1 — A\\ > 2e. Then, there exists 5 > such that 
\\Ap 1- A\\ > 2e + 5. Therefore, there exist <j),ip € Ha, \\4>\\ = Ml = 1, such that 



|(<^, (Ap®l-A)iP)\>2e + -. 
Let Q be a finite-dimensional projection on such that 

11(1- 10Q)(/)|1 < — -, and 11(1- 1(g) Q)^|| < 



8|Uir ^'""-8||A|| 
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Then, 

!((/>, {1 (g, Q){Ap ^ 1)(1 ® - ((^, (1 Q)A{1 Q)V')| > 2e + - - 4^-^P|| . 

2 o\\A\\ 

Since H^dp — Aq\\ < e, this imphes 

\{(P, (Aq g - (1 Q)A{1 <g) > e . 

which contradicts (j3.5|) . 

To conclude the proof, note that for a density matrix in diagonal form, p = Y^f^ pk |Cfc)(?fel S,k, we 
have that id (8) p{A) = J2k Pk^^k ■ Therefore we have 

||n(^)-A|| = II ®ll-^|| < J^pfc||% 1- All <^pfc2e = 2e. 

fc fc fc 

□ 

We now explain a local perturbations perturb locally principle that applies in general to any 
states corresponding to an isolated part of the spectrum of a system of which the dynamics has 
a quasi-locality property expressed by an estimate of Lieb-Robinson type. The basic argument 
can be applied for finite systems or for infinite systems in a suitable representation. For the sake 
of presentation, we consider a systems defined on a metric graph (T,d). To each site x € F, we 
associate a Hilbert space Tlx- For finite A C F, we define 

(3.6) nK = ®n. and Ak = ^B{'Hx) 

xeA 

where B{T-Lx) denotes the bounded linear operators over T-L^- There is a natural way to identify 
^Ao C A\; namely identify each A G with A (g) ]1a\Ao ^ -^A- We can then inductively define 

(3.7) Aloe = U 

Acr 

where the union is taken over all finite subsets of F. The completion of ^loc with respect to the 
operator norm is a C*-algebra, which we will assume to be represented on a Hilbert space and 
assume that a family of Hamiltonians of the form H{s) = H{0) + $(s) on this space satisfies 

Assumption 12. li Additionally, we assume that the Heisenberg dynamics t^^^\ generated by H{s), 
satisfies a Lieb-Robinson bound uniform in s. 

Assumption 3.3. There are constants C{A, B), a > and a Lieb-Robinson velocity f > such that 
for ah s G [0, 1] 

\\[t"^'\A),B]\\ < c(A,5)e-°('^('''PP^'™PP^)-''l*0 
Here, C{A,B) is of a suitable form such as C||yl|| ||i?|| min(| suppyl|, |suppS|). 

Furthermore, we assume that there is a fixed finite subset X CT such that ^'{s) G Ax and 

(3.8) ll^'ll = sup ||^>'(s)|| < oo. 

0<s<l 

The generator D{s) defined in p.Sp and ()2.19p for the local perturbation $(s) is not strictly 
local. However, the fast decay of the weight function W^{t) in combination with Assumption 13.31 
imply that the effect of D(s) is small far away from X. To make this precise, let R > 0, and denote 
by Xr the following 'fattening' of X: 

(3.9) Xr = {x : 3y£X s.t. d{x, y) < R} . 

The following result shows that in the situation described above the unitary U{s) of (j2.3p in 
Proposition 12.41 can be well-approximated by a unitary Vr{s) G Ax^, i-e., with support in Xr. 
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Theorem 3.4 (Local Perturbations Perturb Locally). For any R > 0, there exist unitary operators 
Vr{s) with supp(Vr(s)) C Xr and a constant C, independent of R, such that 

\\U{s)-Vr{s)\\<CG^^\^-:^) 

with G^^^ the subexponential function defined in Lemma \2.6l Consequently, we also have 

(3.10) IIP(I) - VRil)P{0)VR{ir\\ < 2CG^'\^). 

Proof. We begin by defining a local approximation of the self-adjoint generator D(s) starting from 
(j2.19p . Consider the decomposition ^loc = Axr; ® '^txXr ^-iid let Hr = id ^ p for some state p on 
-^r\Xfli s-iid define 

/oo 
dt {t)UR (e"-f^(") (s)e-'*^(^) ) . 
-oo 

Then, for any T > we have the following estimate: 

\\D{s)-Dr{s)\\ < ||$'|| [ dt\W.,{t)\ + \\W4oo rdt||(id-n^)(e'*^Wc^'(s)e-'*^W)||. 

J\t\>T J~T 

For the first term, we apply the bound of Lemma 12.61 part (iv) and for the second term we use 
(j2.13p and Lemma 13.21 and Assumption 13.31 to get 

\\D{s) - Dr{s)\\ < ||$'||2CG(^)(7T) + ^C||$'|| jXle-^^^""^) . 

For the simple choice T = R/(2v), for not too small R, the second term is negligible compared to 
the first , and we obtain 

(3.11) \\D{s) - Dr{s)\\ < C'\\'^'\\G('\^). 
Now, let Vr{s) be solution of 

-^§-/Ri^) = Dr{s)Vr{s) , Vr{0) = 1. 

The claim follows by integrating the estimate (13. lip . □ 

To illustrate this result, we consider the case where the isolated part of the spectrum, Si(s) 
in Assumption 2.1, consists of a non-degenerate ground state energy. Let ipo{s) denote the corre- 
sponding normalized eigenvector and let A G Aw^Xn be an observable supported away from the 
perturbation, whence [j4, VrJ = 0. By applying Theorem 13.41 we immediately obtain 

i(^(.),^^(s))-(^(o),^^(o))i = \{m,u{sr[A,u{s)]m)\ 

= \{m,u{sr[A,u{s) -VR{s)]m)\ 

< 2\\A\\\\U{s)-Vr{s)\\<2G\\A\\G('H:^) 
This estimate clearly expresses the locality of the effect of the perturbation on the state V^(s). 

4. The spectral flow and quasi-locality 

The main goal of this section is to prove that the spectral flow deflned in terms of the unitary 
operators U{s), as in Proposition 12.41 satisfies a Lieb- Robinson bound. This is the content of 
Theorem 14.51 below. In Section [4.11 we introduce the basic models to which our result applies and 
state Theorem 14.51 Our proof of Theorem 14.51 demonstrates that the claimed estimate follows from 
a Lieb-Robinson bound for time-dependent interactions. We state and prove a general result of this 
type, see Theorem 14.61 in Section 14.21 The remainder of Section |4] is used prove that Theorem 14.61 
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is applicable in the context of the spectral flow. Section 14.31 contains a technical lemma, and 
Section 14.41 finishes the proof. 

4.1. The set-up and a statement of the main result. The arguments we provide in Section U] 
apply to a large class of models. In this subsection, we describe in detail the assumptions necessary 
to prove a Lieb-Robinson bound for the spectral flow. 

We will consider models defined on a countable set T equipped with a metric d. Typically, F 
will be infinite, e.g., T = W. In the case that F is infinite, we require some assumptions on the 
structure of F as a set. First, we will assume a uniform bound on the rate at which balls grow, i.e., 
we assume there exist numbers k > and v > for which 

(4.1) sup \Brix)\ < Kr" , 

where is the cardinality of the ball centered at x of radius r. In addition, we will assume 

that F has some 'integrable' underlying structure. We express this property in terms of a non- 
increasing, real-valued function F : [0, oo) (0, oo) that satisfies 

i) uniform integrablility: i.e. 

(4.2) ||F|| = sup VF((i(x,2/)) < oo 
and 

ii) a convolution condition: i.e., there exists a number Cp > such that given any pair x,y 

(4.3) Yl ^))^(d{z, y) < CFF{d{x, y)) . 

zer 

For the case of F = W , one possible choice of F is given by F{r) = {l + r)~'^'^~^^\ The corresponding 
convolution constant may be taken as Cp = 2^^^^ Sxer -^(kl)- 

Lastly, we need an assumption on the rate at which F goes to zero. It is convenient to express 
this in terms of the sub-exponential function Ua introduced in Lemma 12.51 We suppose that there 
exists a number < (5 < 2/7 such that 

(4.4, s?i|y<~- 

Clearly, if F = and F(r) = (1 r)-(''+i), then g3D holds for every < 5 < 2/7. 

The following observations will be useful. Let F : [0, oo) (0, oo) be a non-increasing function 
satisfying (|4.2p and (|4.3p . For each a > 0, the function Fa{r) = e~°'^F{r) also satisfies the properties 
(j4.2p and (j4.3p with \\Fa\\ < and ^ Cp- In fact, more generally, if g is positive, non- 
increasing, and logarithmically super-additive, i.e., g{x + y) > g{x)g{y), then Fg{r) = g{r)F{r) 
satisfies (|4.2p and (|4.3p with \\Fg\\ < g{0)\\F\\ and Cpg < Cp- For brevity we will write Fa to denote 
the case g{r) = e""*". Other functions g will be used later. 

Recall the general quantum systems corresponding to F on which our models will be defined. 
As in Section [3l we associate a Hilbert space T-La and an algebra of observables Aa to each finite 
set A C F, see (13. 6p . and similarly define Aioc as in ()3.7p . In this case, the models we consider 
are comprised of two types of terms. First, we fix a collection of Hamiltonians, which we label by 
(-ffA(O))^, with the property that for each finite A C F, H\{0) is a densely defined, self-adjoint 
operator on 7i\. Next, we consider a family of interactions ^{s) parametrized by a real number 
s. For each s, the interaction ^(s) on F is a mapping from the set of finite subsets of F into ^loc 
with the property that ^{X,s)* = ^{X,s) S Ax for all finite X CT. It is convenient to write 
$(X, s) = ^x{s)- A model then consists of a choice of {Ha{0))^ and a family of interactions <I>(s) 
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over r. Given a model, we associate local Hamiltonians to each finite set A C F by setting 

(4.5) Ha{s) = Ha{0) + J2 ^xi^) 

where the sum is taken over all subsets X C A. For notational consistency, we will assume that 
^x(O) = for all X. With s fixed, the sum in (j4.5p above is finite for each such A C F, and thus 
self-adjointness guarantees the existence of the Heisenberg dynamics, i.e., 

(4.6) 7-f^^'^(^) = e^*-^A(s)^g-itHA(s) foj, ^ £ and t E M, 

which, again for fixed s, is a one-parameter group of automorphisms on ^a- 

To prove the results in this section, we need a boundedness assumption on the family of inter- 
actions. We make this precise by introducing a norm on the interactions $(s) over F, with respect 
to any non-increasing, positive function F satisfying (j4.2p and ()4.3p . as follows: 



(4.7) \Mf= sup — r- sup||«>z(s)|| <oo. 

x,y£Z 

The sum above is over all finite sets Z CT containing x and y, and we will often abbreviate || • \\f^ 
by II • lla- On occasion, we will use ||<I>(s)||i? for the norm of $(s) at fixed s, i.e., the norm defined by 
dropping the supremum over s in (|4.7p . The following lemma states some simple bounds in terms 
of \\^\\f that we will frequently use. 

Lemma 4.1. Let ^{s) be a family of interactions over F for which \\^\\f < oo for some non- 
increasing, positive function F satisfying and |^.3| ) above. Then, for any finite A C T , we 
have 



\F 



(4.8) 11^^(^)11 ^ ^Wll*^! 

XCA: 
x&X 

(4.9) Yl \\^x{s)\\ < F(0)||$||f|A|. 



XcA 



Proof. For x G F we have 



Ej|.x(.)ii<™pm.,.))|:i«£no)ii*ii. 

x£X x,y&X 

where we have used the definition of the norm (|4.7p and the monotonicity of F. Using this estimate, 
for any finite subset A C F, we obtain the bound 



Y \\'^x{s)\\ < E E ^ ^(o)ii^iif|ai 



XcA xeAXcA: 

xex 



□ 



We will also require the interactions to be smooth with bounded derivatives. More concretely, 
let $(s) be a family of interactions over F for which, given any finite X cT, ^x{s) is differentiable 
with respect to s. In this case, we define a corresponding family of interaction d^{s) over F by the 
the formula 

d^x{s) = \X\^xis) for each finite X C F. 
We now state the main assumptions of this section. 
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Assumption 4.2. We will assume that the interactions $(s) are differentiable with respect to s. 
More specifically, we assume that for each finite X C T, G Ax for all s. In addition, we 

suppose a uniform estimate on the norms of these derivatives as s varies in compact sets. For 
concreteness, we will assume that the domain of s- values is [0, 1], and suppose that there exists a 
number a > for which 

\\d^\\a < OO. 

Assumption 4.3. We will assume that for every finite A C T, the local Hamiltonian H\{s) has 
a spectrum that is uniformly gapped. More precisely, the spectrum of H\{s), which we will de- 
note by S('^)(s), can be decomposed into two non-empty sets: S(^)(s) = T.[^\s) U J^i^\s) with 
d{T,[^\s),T,^\s)) > 7 > 0. In particular, the positive number 7 is independent of s € [0, 1] and 
finite A C r. We also suppose that there exist intervals I{s), with endpoints depending smoothly 
on s, for which C I{s). 

In typical applications, the set will consist of the ground state and (possibly) other 

low- lying energies, but this is not necessary. 

Given Assumptions 1412) and 14131 the results of Section [2] apply to the local Hamiltonians H\{s). 
We need a further assumption in order to state the main result of this section. 

Assumption 4.4. We will assume a uniform, exponential Lieb-Robinson bound. In fact, we assume 
that there exists an a > and numbers Ka and Va such that 

(4.10) ||[rf'^(^)(A),i?]|| <K„P||||S||e'^^'"l*l Fa{d{x,y)) 

xex,yeY 

holds for all A G Ax, B € Ay, and t G M. Here, as above, Fa{r) = e~'^^F{r), and we stress that 
the numbers Ka and Va are each independent of both A and s. 

Estimates of the form (j4.10p have been demonstrated for a number of models, see e.g. [H], and 
references therein, for a recent review. Here we assume it holds for a class of models, and as a 
consequence, we get Theorem 14.51 below. 

As indicated above, given Assumptions 14.21 and 14.31 the results of Proposition 12.41 applv to H\{s) 
for each finite A C T and s € [0, 1]. In this case, there are unitaries U\{s) in terms of which we 
define the following spectral flow: 

(4.11) a^{A) = Ua{s)*AUa{s) for all A € and < s < 1 . 

The main result of this section is a Lieb-Robinson bound for the spectral flow, which is formulated 
with the aid of a function Kjr defined as follows: 

(4.12) F^{r)=uJ^r]F(^r 
where 

(4.13) ujx) - . , . . 

^ ' Uf^{x) otherwise. 

Since F is uniformly integrable over T and u^{r) < 1, Fx^/ satisfies ()4.2p . Moreover, Fq, also satisfies 
(j4.3p . In fact, it is easy to check that is positive, non-increasing, and logarithmically super- 
additive. The Lieb-Robinson velocity in the following theorem also involves the norm H^'Hf^, of an 
interaction ^ defined later in this section (see (|4.46p ). 



8Va J \8Va 



n^(e^) for < r < e^, 



Theorem 4.5. Let Assumptions \4-^ \4-3[ and \4-4\ hold. Then, 
(4.14) ||[a^(A),B]|| < 2||A||||B||min 



l,g{s) F^(.d{x,y)) 

xeX,yeY 
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for any A € Ax, B S Ay, and < s < 1 and g is given by 

e2|l*llF^Cf^|t| _ 1 ifd{X,Y)>0, 



(4.15) CF^-g{t) 



e 



'iWnF^CF^ltl otherwise. 



The number Cp^ is as in |^.3[ ) and our estimate on ||^||_f^ is discussed in the next subsections. 

4.2. Lieb-Robinson bounds for time-dependent interactions. The estimate ()4.14p in the 
statement of Thereom 14.51 can be understood as a Lieb-Robinson bound for the spectral flow. In 
this section, we demonstrate that Lieb-Robinson bounds hold for a large class of time-dependent 
interactions. As in the previous section, we assume that our models are defined on a countable set 
r equipped with a metric. Let $t denote a family of interactions over F, and, for convenience, we 
win assume that t G [0, 1]. Thus, for every finite X CT and each t G [0, 1], ^t{X)* = ^t{X) G Ax, 
and we will often write ^t{X) = ^x{t). 

In this case, corresponding to each finite A C L, there is a time-dependent local Hamiltonian 
which we denote by 

(4.16) HA{t) = ^x{t) . 

XcA 

We will assume that, for each finite A C F, HA{t) is a strongly continuous map from [0, 1] into Aa- 
In this case, see e.g. Theorem X.69 [U], it is well-known that there exists a two-parameter family 
of unitary propagators UA{t,s) with 

(4.17) j^UA{t,s) = -iHAmAit,s) and C/a(s,s) = ]1, 

the above equation holding in the strong sense. The Heisenberg dynamics corresponding to HA{t) 
is then defined by setting 

(4.18) rl'iA) = UA{t, OyAUAit, 0) for all A e Aa ■ 
The following Lieb-Robinson bound holds. 



Theorem 4.6. Let F be a non-increasing, positive function satisfying |^.g[ ) and ( |^.3| ) and suppose 
that the interactions satisfy 

(4.19) \Mf= sup / Yl ^^'P < 

x,yar F{d{x,y)) ^ o<t<i 

x,yGZ 

Then, for any subsets X, y C F, ^ G Ax and B G Ay the estimate 

(4.20) ||[r/^(A),fi]|| < 2PIIIIS 

where the function g may be taken as 

(4.21) CF-g{t) = 



mm 



l,g{t) Y F'id{x,y)) 



e2||*||FC^|i| _ 1 ifd{X,Y)>0, 
^2\\<s>\\FCF\t\ otherwise, 

and the number Cp is as in ( |.^.3[ ). 

Proof. Let X,Y C F be finite sets. Take A C F finite with X \JY C A. Define the function 
/ : [0, 1] ^ >iA by setting 

(4.22) f{t) = [UA{t,^TUx{mAUx{t,QTUA{t,Q),B] = [t^ {f^{A)),B] , 
where we have introduced the notation {A) = Ux{t,Q)AUx{t,Q)* ■ Denoting by 

(4.23) = {Z C A : znx / 0,Zn 0}, 
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the surface of X, a short calculation shows that 

fit) = i[Ti^{[HA{t)-Hx{t),fi'{A)])B] 



ZCA: ZeA: 

zes^ Z£S^ 



As the first term above is norm-preserving, see e.g. j36| . the inequality 

(4.24) ||[r,^(f/^(^)),i?]||<||[Ai?]|| + 2P|| Y r\\[rH^z{s)),B]\\ds 

Jo 

follows. Consider now the quantity 



(4.25) C5(.Y,.)= .up 

AeAx- 



\A\ 



It is easy to see that 
(4.26) C^(X,t) <C^(X,0) + 2 V sup ||$z(r)|| / C^iZ,s)ds. 



From here, the argument proceeds as in the proof of Theorem 2.1 in [36]. □ 

4.3. Some notation and a lemma. In this subsection, we prove a technical estimate needed in 
our proof of Theorem 14.51 The objective is to show that the s-dependent generator of the unitary 
flow U\{s) has the structure of a bonafide short-range interaction. In Theorem 13.41 we showed that 
each term of the perturbation, i.e., $x(s) for a given X, leads to a term in the generator that can 
be well approximated by local self-adjoint operator supported in Xji with almost exponentially fast 
decay of the error as a function of R. A projection IIx^ : ^a •^x^ was used to accomplish 
this. In this subsection and the next we apply the same procedure to show that the differences 
between successive approximations can be summed leading to a decomposition of each term in the 
generator as a telescopic sum of finitely supported terms. To define the terms in this decomposition 
we need a family of projection mappings {Ilx)xcA^ ^^^'^ decomposition we obtain will depend 
on the choice of this family. It will be convenient to choose a family which is compatible with the 
embeddings ^a^, C Aa, for Aq C A, and such that each of the Tlx are continuous in the norm 
and weak topologies on Aa- We will therefore choose a family of normal states on B{'Hx), or 
equivalently, a family of density matrices, {px)x£r define a product state on Ax'= by 

setting px': = 'S>x&r\x P^- Then, for any finite X C A, we define 

(4.27) = idAx ^ Px-\aa ■ 

Here, id^^^. is the identity map on Ax- can be considered as a map — ^ Aa with ranllx C Ax- 
We let the dependence of Tlx on the px be implicit. All our estimates will be uniform in the px- 

Similarly, the interaction ^a(s) we define in the next subsection depends on the choice of px, but 

the estimates on its decay will not, and the unitary flow generated by these interactions also does 

not depend on the px- 

Fix a finite set A C F. For any X CT and n > 0, denote by 



(4.28) 



Xn = {zer: d{z,X) < n} , 
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J —oo ^ ' 



where d{z,X) = m.mx£x d{z , x) . Keeping with the notation from the previous subsection, for any 
A G Ax we set 

(4.29) 
and 

nx„ (rf-^'\A)) W,(t) dt - / nx„_, {r!''^'\A)) W,{t) dt 

-OO "/ — oo 

for any n > 1. Since A is finite, A^(^, s) = for large n. Moreover, it is clear that supp {A'^{A, s)) C 
Xn n A. In our proof of Theorem 14.51 we will use that 



(4.31) 



/oo °° 
T^'^'\A)W,it)dt = Y,^l{As) 



where the series is actually a finite sum. In fact, the following estimate is also important. 

Lemma 4.7. Under Assumptions \4-^ and \4-4l let A CT be a finite set. For any X C A, A & Ax, 

and integer n > 0, 

(4.32) ||AX(^,s)|| <2p||min[||VF^||i,|X|G(n-l)] 

where 



(4.33) 



and is as in Lemma \K 



G{n) = 4/, ( ^\ + :^e— /2 



Proof. It is easy to see that 

(4.34) ||A^(^,s)|| < Pll llW^lli and ||AX(A, s)|| < 2P|I . 

A better estimate in n is achieved by inserting and removing an identity. In fact, we need only 
estimate the norm of 

(4.35) 



/OO , , 

(nx„-id) [r^'^^'\A))w,{t)dt. 
-oo 



To do so, we follow the same strategy as in the proof of Theorem 13.41 By Assumption 14 . 4 ^ we know 
that 



(4.36) 



{A),B 



< A'J|F|||X||M||e"'^''l*le-""||5| 



for all B € Ax'^ ■ Hence, for any T > 0, we have that 



(n 



l\t\<T 

(4.37) 

using Lemma 13.2^ whereas 



id)(rf^(^)(A))l^,(t)dt 



< 



I / (nx„-id)(rf^(^)(^) 

^ J\t\<T ^ 



dt 



< Ka\\F\\\X\\\A\\ 



dt e 



(4.38) 



/ (nx„-id)(rf^(^)(^))i^,(t)dt 

J\t\>T ^ ' 



< A\\A\\ly{T) . 



The choice of T = n/2va yields an estimate of the form 



(4.39) 



r {Ilx^-id)(r!'^^'\A))w,{t)dt <A\\A\\lJ^)+^^\X\\\A\\e---/\ 



The bound (|4.32p readily follows. 



□ 
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As indicated by the proof above, a stronger inequality is true. We have actually shown that for 
every n > 1, 

(4.40) \\Al{A,s)\\ < 2\\A\\mm[\\W4i,GA{n - 1) + GA{n)] 
where 

(4.41) G,(„, = 2/,(^)+M3|x|.-W.. 

For the arguments we use below, it is convenient to extract a decaying quantity that is independent 
of the given observable A and use the monotonicity of G. This explains the form of the bound 
()4.32p appearing in Lemma 14.71 

4.4. The proof of Theorem 14.51 In this subsection, we prove Theorem 14.51 The basic idea is 
that Theorem 14.51 follows from a Lieb-Robinson bound for time-dependent interactions, see e.g. 
Theorem 14.61 in Section 14.21 To see that such a result is applicable, we demonstrate that the 
generator of the spectral flow can be written as a sum of local interaction terms which satisfy an 
appropriate decay assumption. This is the content of Theorem 14.81 below. 

Under Assumptions 14. 2l and 14. 3l we have defined (for each finite A C F) a spectral flow by setting 

(4.42) a^{A) = Ua{s)*AUa{s) for all A G Aa- 

In fact, the unitary Ua{s) is the one constructed in Proposition 12.41 and as a consequence of 
Corollary 12.81 we know that C/a('S) is generated by 

/oo 
rf^(^) {H'^{s))W,{t)dt 
-oo 

/oo 
rf^(^) {<^',is))W,it)dt. 

Here 7 is as in Assumption 14.31 and appears in Lemma 12.61 . The previous subsection demon- 
strated that each term 

/oo 
rf (^'z(^)) W,{t) dt = Y, ^m'zis), s) 

n=0 

where the series is actually a finite sum. Combining ()4.43p and ()4.44p above, we write 

00 

(4.45) Da{s) =Y.11 ^li^'zi^)^') = E ^a(^' ^) ' 

ZcAn=0 ZcA 

where 

(4.46) ^a{Z,s) = Y^ ^li^Y(s),s). 



n>0 YcA: 

Yn = Z 



It is important here to note that supp(^A(.2^) s)) C Z, i.e., the s-dependent, interaction terms 
^AiZ,s) are strictly local. The following estimate holds. 



Theorem 4.8. Let Assumptions \4-^ \4-S\ and \4-4\ hold. Then, there exists a function Fq, satisfying 
\4-^ fl^t^ |-^-3| ) such that 

(4.47) II^aIIf* = sup — -J V sup \\^a{Z,s)\\<oo. 

x,y£Z 

Here we note that the function is independent of A. 

It is now clear that Theorem 14.51 follows from Theorem 14.81 via an application of Theorem 
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Proof. In the argument below, it is convenient to set a > to be the minimum of the a's whose 
existences are guaranteed by Assumptions 14.21 and 14.41 
We begin by re- writing the quantity of interest. Clearly, 



(4.48) 



sup ||^a(^,s)||< Y1 ""^P II^A(^y(s),s 



0<s<l 



y,7i>0: 



0<s<l 



and so 



ZCA: 



0<s<l 



(4.49) 



sup ||^a(Z,.)|| < 5] 2; sup \\Al{^'y{s),s)\\ 

x,y(^ZY„=Z 

= 5^ 5Zlnd[x,y sup \\AI{^'y{s), s) 

YcAn>0 ^^'^^ 

= E E 



ycA: n>0 



0<s<l 



+ J2 E lnd[{x,y}ny„';_i^0] sup \\Al{^'y{s),s 



m=l YCA: 



n>m 



0<s<l 



The first equality above follows from the observation that 



ZCA: Y,n>0: 

x,y&ZY„=Z 



Z<ZhY(ZKn>0 



= El"d[y„ = Z] Ind [x,yGy„] 

YcAn>0 UcA 

(4.50) = ^ ^Ind[x,y , 

YcAn>0 

while the second is a consequence of the fact that for any pair x, y 

(4.51) E = E + E E [{x, y} n y^_i / 0] . 

ycA rcA: m>l rcA: 
x,y&Y X,y&rn 

The first sum on the right-hand-side of (j4.49p is easy to bound. In fact, using Lemma 14.7^ it is 
clear that 



(4.52) 



sup \\AI{^'y{s),s)\\<2\Y\ sup ||$y(s)||G(n- 1) 

0<s<l 0<s<l 



where G is as in ()4.33p with G(— 1) set to be HW^Hi- Thus, 

Y E ^^P l|AX(^'y(^),^)ll < 2YG{n-l) ^ \Y\ ^^P ll^y(^) 



YCA: n>0 
x,y(iY 



0<s<l 



n>0 



YCA: 
x,y(iY 



0<s<l 



(4.53) 



< 



2\\dnaFa{d{x,y))YG{n-l). 



n>0 



From the estimates in Lemma 12.61 it is clear that G is summable. 

For the remaining terms in (j4.49p . we use the following over-counting estimate: 



(4.54) 



Y Ind[{x,y}ny,^,_i/0] < E E E 



YCA: 

x,y&m 



yi&B^{x) y2&B^{y) YcA: 

yi,y2&Y 
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Combining (f02]) with (fi3i]) . we find that 



^ Ind[{x,y}ny^_i/0] sup \\Ali<P'y{s), s)\\ 

\i . \ 0<s<l 

m>l 1 cA: n>m — — 

E E E m^supji$v(^)iiE^(^-i) 

m>ly^eBm{x)y2&Bm(y) ^CA: n>m 

(4.55) <2\\dnaY.G{m) E E ^-(^i'^^)), 

m>l yi&Bmix) y2&B,n{y) 

where we have set 

(4.56) dim) = Y^Gin-l). 

n>m 

We now perform a rough optimization over m > 1. Take < e < 1 and declare uiq = mo(e) > 
to be the largest integer less than (1 — e)d{x, y)/2. We claim that, for m < niQ and yi and y2 as in 
(j4.55p above, ed{x,y) < d{yi,y2)- This follows from 

(4.57) d{x, y) < d{x, yi) + d{yi,y2) + d{y2,y) < d{yi,y2) + 2m< d(yi, ^2) + 2mo , 
and the choice of rriQ. In this case we have 

mo+l mo+1 

E ^("^) E E P-id{yi,y2)) < G{l)Fa{ed{x,y)) E \B„^{x)\\B„,{y)\ 

m=l yi€B,n{x) y2&B^{y) m=l 

mo+1 

(4.58) < K2G(l)F,(ed(x,y)) E 



m=l 



where we have used (j4.ip . 

The remaining terms we bound as follows. 

E ^("^) E E F-(d{yi,y2)) < \\Fa\\ E \Bm.ix)\G{m) 

m>mo+l yi&B,n{x) y2&B,n{y) m>mo+l 

(4.59) < KllFall E "^''G'(m). 

m>mo+l 

Now, from the definition of G, 



(4.60) E E ™" E „„„ 

m>mo + l m=mo+2 n=m—l ^ ^ ' 

and the sum 



g-an/2 



00 00 

u —am/2 



(4.61) E E e-""/2^ 6-^/2 ^e-'^^/2. ^ ^ 

m=mo+2 n=m— 1 J/>0 m=»no+2 
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decays exponentially in mo- Using the results in Lemma 12.51 and 12.61 we find that 

OO OO / ^ „00 OO/^10 / 

E E ^. G^^ f E E £ £ 

m=mo+2 n=m-l \ "/ ' m-mo+2 n=m-l ^ \ " 



< 



< 



< 



2VaC 



7^ 



m=mo+2 \ ^ \ 

' ' + 1 y u2/-j{y)dy 



^3 / 7(?I1Q+1) \ 



This proves that 



V —am/2 



sup ||^'a(^,s)|| < CiF,(ed(x,y))(mo + l)''^+^ + C2 m"e 

zcA: '^^''^1 m=mo+2 

for some number p depending only on v. Since 2mo < (1 — e)d{x,y), it is clear that the final term 
above decays the slowest in d{x,y). Thus we have shown that 

(4.64) „^^P^ \\^AiZ,s)\\ < [^d(x,y) + l^J u,/r , 



ZCA 
x,y£Z 



for each < e < 1. For concreteness, take e = 1/2. With 5 > as in ()4.4p and any < 5' < 2/7 — 6, 
we will set fi = 2/7 — 6 — S' > and see that 

Y sup \\^a{Z,s)\\ < C'(l^d{x,y)yu,/J^d{x,y)] 

0<^*<l \8Va J ' \8Va ) 



ZCA: 



(4.65) < C'"u^\^d{x,y)y\^d{x,y)y 

With the definition of given in (j4.12p , this completes the proof of (|4.47p . □ 

5. Existence of the thermodynamic limit and gapped quantum phases 

The Lieb- Robinson bound for the flow given in Theorem 14.51 of the previous section, can be 
used to obtain the thermodynamic limit of this flow defined as a strongly continuous cocycle of 
automorphisms of the C* -algebra of quasi-local observables. The standard setting is the same as 
in the previous section, but we now assume that the Hilbert spaces 'H.x associated to each x € T, 
are all finite-dimensional. The C*-algebra of quasi-local observables Ay is then obtained as the 
completion with respect to the operator norm of A\oc- 



(5.1) = Axoc = U 



ACF 
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If Hx is allowed to be infinite-dimensional it is typically necessary to work in the GNS representation 
of a reference state in order to have a well-defined thermodynamic limit. Such an approach was 
used in [38] to define the dynamics of an infinite lattice of anharmonic oscillators. In order to avoid 
the need for additional technical assumptions, for the remainder of this section we restrict ourselves 
to quantum spin systems, i.e., the case of finite-dimensional 1-ix- It is not necessary, however, that 
dYml-ix is independent of x or even uniformly bounded. 

This section has two subsections. In the first, we prove that the finite volume spectral fiows, 
defined as in (j4.1ip . have a well-defined thermodynamic limit. With these results in hand we can 
then, in the second subsection, complete the proof that gapped ground states connected by a curve 
of quasi-local interactions satisfying a suitable norm condition are equivalent under a quasi-local 
automorphism, in finite volume as well as in the thermodynamic limit. But first we describe in 
detail the class of systems to which our main result applies. 

The systems under consideration here have finite dimensional local Hilbert spaces. In this case, 
we can make a convenient choice of the projection map introduced in Section 14.31 and needed for 
the application of Lemma 13.21 namely the natural extension of the partial trace. For any finite 
subset A C r, we define the conditional expectation IIa : — > -^A as 

IIa = id^A ® ^Ak'^ , 

where for A' C T, 

TA,, = (g) r^. , = d-^Tr^^ 

a;eA' 

is the normalized trace over A\i. In particular, for any Z C A C T, the subprojections 

nA,z = nz|_4^ 

form a consistent family, namely for any A € Ax, with Z^X d Am C A„ C F, they satisfy 

(5.2) nA„,z(^) = ^A^M^) 

and the first index may be dropped. 

Let F be a countable set equipped with a metric and a function F satisfying (j4.2p and (j4.3p . For 
s € [0, 1], let ^{s) be a family of interactions, differentiable in s, for which there exists a number 
a > so that 

(5.3) \Ma+\m\a<00. 

where the norm is defined in the paragraph containing (j4.7p . 

Our proof of the existence of the thermodynamic limit requires some assumptions on the sequence 
of finite volumes (A„)^ on which the spectral flows are defined. Let (A„)^ be an increasing sequence 
of finite sets which exhaust F as n — )■ co. For convenience, we will regard the parameter n as 
continuous with the understanding that, for any n > 0, A„ = A[„], where [n] denotes the integer 
part of n. We will assume that there exist positive numbers bi, b2, and p such that 

(5.4) (i(Am, A.^) > bi{n - m), and |A„| < b2nP . 

We assume that there are finite intervals I{s), smoothly depending on s € [0, 1] such that, for all 
n, the finite-volume Hamiltonians ifA„(s) = X]^cA„ ^{^^^) have one or more eigenvalues in I{s), 
and no eigenvalues outside /(s) within a distance 7 > of it. 

Let us summarize the results of the previous sections, given these assumptions. If PA„is) 
denotes the spectral projections of H\^{s) on I{s), then there is a cocycle a^" , the dual of 
which maps -Pa„(0) to Pa„{s) for all s € [0,1]. Its generator has a local structure given by 
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-^A„('5) = X^zcA„ ^A„(-2^>s) where the interactions ^'a„(s) decay almost exponentially in the fol- 
lowing sense, 



(5.5) 



|^a„||f*= sup J" 



X] ll*A„(^,s)|| < oo. 



ZcA, 



0<s<l 



uniformly in n, where satisfies again the uniform integrability and convolution property for a 
constant C^. Our estimates in Section S] demonstrate that a possible choice of is given by ()4.12p 
which decays sub-exponentially. 

5.1. Thermodynamic limit for the spectral flow. In order to prove the existence of the ther- 
modynamic limit of the spectral flow a^, it is convenient to recall an estimate from the proof of 

, as proven e.g. in [36]. 



the existence of the thermodynamic limit of Heisenberg evolutions 
In fact, assuming that ||$||a < oo, the following bound is valid. 

Take finite sets X C C A„,. Note that for any A € Ax, each s € [0, 1], and any t G M, 



Ma(s) 



(A) 



(5.6) 



< 



< 



< 



E 

ZCA; 

Ka\\A\ 



ZCA„: 
ZnAn\Am¥=<l> 



\t\ 



''^zis),Ti;^t-^'\A) 



dr 



av„ 



ava\t\ 



1) E 11^^(^)11 E PMi^^^)) 



ZCAn: 

ZnA„\A„^0 



^'^"^"C^J|c^|U(e-"l*l-l) E Fa{d{x,y)). 



av„ 



yeA„\A„ 



Since Fa is uniformly integrable, this proves that the sequence (^t^'^"^^\a)^ is Cauchy. We will 



denote the limit by ' (A), and observe that it satisfies 
(5.7) 



rp(A)-rf^-(^)(^) 



ava 



1) E F'aidix,y)), 



xex 
yeT\A„ 



uniformly for s € [0, 1]. 

The following analogue of Lemma 14.71 will be useful. Recall the definitions of A'^(A,s) from 
(j4.29p and (|4.30p . Define similarly Ar{A,s) with t^'^{A) replacing t^'^^^^A) as appropriate. 

Lemma 5.1. Let A C L be a finite set. For any X C A and A € Ax, 
(5.8) \\Al{A,s) - A'^{A,s)\\ < 4||A||min |A:| VG(n - l)K{d{X, A^)) 

where G is as in 1^4-33 ) of Lemma \4- ?| and 



(5.9) 



K(x) 



AF 



X 

2^ 



2, ,2 



Proof. A uniform estimate, as shown in Lemma [4.7| clearly holds for n = 0. We need only consider 
n > 1. Using the consistency of the mappings IIx^, the difference A^{A,s) — Ap(y4, s) can be 
written as a difference of two terms. As such, we need only bound the norm of 



(5.10) 



/OO , . 

(nx„ - id) U'^A) - tI^'^'Ha)) w,{t) dt , 
-oo 



24 



S. BACHMANN, S. MICHALAKIS, B. NACHTERGAELE, AND R. SIMS 



By Assumption 14.41 t^^^^'^ satisfies a Lieb-Robinson bound uniform in A and s. In this case, the 
Hmit Tf does as well. Arguing then as in Lemma |4.7| it is clear that 



r (n^„ - id) (r^'^iA) - Tf^^'\A)) W,it) dt 



(5.11) 



< 2|X|||yl|l 4/. 



'7 



n 

2^ 



+ 



Ka\\F\\ 



-an/2 



Since the projections Tix„ are norm one maps, we may also argue using the thermodynamic 
estimate ()5.7p . In fact, 

(nx„ - id) U'^iA) - r^'^'\A)) W,{t) dt <2 rf'^iA) - Tf'^'\A) \W,{t)\ dt . 

Now for |i| < T, we have that 

2 I Tl^\A)-r^'^^'\A) \W,{t)\dt < ^^^!^I^\\A\\ F,{d{x,y)) f e-^l^'dt 

J\t\<T aVa ^ J\t\<T 



'\t\<- 

(5.12) 

whereas for It I > T, the bound 



yer\A 



< 2i^aq^J'^IU p|||^|||j,||g-ad(X,A^)^a..T^ 



(5.13) 



2 / rp(A) - t^^^'\A) \W^{t)\ dt < 8P|l/^(r) , 

^|t|>T 



is clearly true. In this case, the choice T = d{X, K^) / [2va) yields the estimate 



(n 



id)[Tl'\A)-Tr'"{A))W,{t)dt 



(5.14) 



< 2IXIMII AL 



fd{X, A^) 



'7 



2v„ 



+ 



-ad{XA'')/2 



Combining the results from ()5.1ip and (|5.14p . as well as the bound corresponding to IIx^.i, the 
estimate ()4.32p follows. □ 



We can now state and prove the existence of the thermodynamic limit for the spectral flow a 
Recall that for any finite sets Z C A C T, we have defined 



(5.15) 

By analogy, set 
(5.16) 



^j,{Z,s)= AX(a>V(s),s). 



Y,n>0: 
Yr,=Z 



^r{Z,s)= J2 AF($'y (.),.) 



y,Ti>0: 
Yr, = Z 



We will show later in this subsection that the s-dependent interaction ^'r(s) is the limit as A F 
of ^'a(s). First, we show the existence of the limiting spectral flow in Theorem 15. 2[ Then, we 
argue that it is also the limit of the automorphisms generated by finite volume restrictions of the 
limiting interaction ^'r(s). 

Theorem 5.2. Let (o^")^ denote the sequence of flows associated with the sets A„ C F. Then 
there exists a flow defined on the quasi-local algebra Ar such that for all A G A\oc, 



hm \\at-iA)-al{A)\\=0., 
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uniformly for all s £ [0, 1]. 

Proof. We begin by noting that the strong Umit of an automorphism is automatically an auto- 
morphism and that convergence of a sequence of automorphisms cj„ — )• a, is equivalent to the 
convergence of the inverses to the inverse automorphism, i.e., a^^ a~^. Using these observations 
and by standard completeness arguments it is therefore sufficient to establish that for all A E ^loc; 
the sequence {a^")~^{A) is Cauchy. Without loss of generality, we assume that A e Aaq and we 
use the notation a^" = (a^")~^. Then, for n > m, define 

/(.)=d^(^)-d^(^). 

and observe that 

fis) = i[D^„{s),a^-{A)]-i[Dj,^is),dj-iA)] 

= z[DA„{s),f{s)] + ilDA^s) - DA„(s),a^(A)] . 

Hence, 

(5.17) = ||/(.)|| < r\\[DA,Sr)-DA^{r),a^-{A)]\\dr. 

Jo 

We will show that the right-hand-side goes to zero as n, m ^ oo. 
We begin by writing the difference as 

DA^r) - DA^ir) = Yl ^A.(^'0 + E (^A„(^,r) - ^A„.(^,r)) . 

ZCA„: ZcAm 
Zn{A„\A„)^0 

For the first term, the Lieb- Robinson bound of Theorem 14. 5 1 which clearly applies to a^'" as well, 
yields 

mA^{Z,r),di^-{A)]\\<2\\A\\\\^AAZ,r)\\g{r) J] F^{d{x,y)) . 
After summing over Z and integrating, we find that 



r Y mAAZ,r),d^-{A)]\\dr 
zn(A,AA™.)7^0 

<2\\A\\ r Yl \\^AAZ,r)\\gir)dr ^ ^^idix,y)) 

° ZcAn. x£Ao,y£Z 

zn{A„\A„)^0 

<2\\A\\ r gir)dr ^ E sup ||^aJZ, r)|| ^ ^^K^' y)) 
° yeA„,^eA„\A„zcA„:°-'''-^ xeAo 

z,y(^Z 

<2\\A\\\mC<, r g{r)dr ^ E^^^^^^'^)) 

° 2eA„\A„i xeAo 

which vanishes as m < n ^ oo by the uniform integrability of . 

To control the second term, we arrange the set of subsets of A^, which we denote by V{Am), as 
a union of three sets: V{Am) = "Pi U 'P2 U "Ps where 

(5.18) Vl = {ZG V{kra) : Z C A^/3}, V2 = {Z £ V{km) : Z C k^m/z}, 
and 

(5.19) P3 = e V{Km) ■■ Z n A„/3 / and Z n A^^^/3 / 0}. 
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We first sum over Vi- Repeating the argument we used above, in particular using the uniform 
Lieb- Robinson estimate for both ^\^{Z,r) and ^\^{Z,r), we find that 

r E II [^A„(^,r)-^A^(Z,r),d^™(^)] II < 4||A|| ||*||C^ r g{r) dr J] J] F^{d{x,z)), 
•^0 •''0 ^gAc^^^^gAo 



zeVi 

and this bound decays to zero as m — )■ oo. 

We next estimate the sum over 7^2- We begin by trivially bounding 



< 2\\A\ 



(^A„(Z,r)-^r(^,r)) 



(5.20) 



+211^1 



(M/r(Z,r)-M/A„,(^,r)) 



where we are using the notation from ()5.16p . Each of the terms on the right-hand-side above will 
be estimated similarly. In fact, note that 

(5.21) 5] (M/r(Z,r)-^A„(Z,r))= ^ J] ^ (AF($V(0, " AA„(<I>V(r), r)) 



implies a bound of the form 



iMZ,r)-^A^iZ,r)) 

Z&V2 



^CA2„/3 n>0 ycT: 

Yn — Z 



< Y H ||AB($V(r),r) -Aa„ ($'y(r),r) 



n>0 YCF: 



(5.22) 



< AYVG{n-l) Y E 1^1 ll'^VWIIV^(^^(W'^ 

< AWd^aFaiO) Y VG{n - 1) • \A2m/3\VK{bim/3) . 



n>0 



Since |A2m/3l ^ b2{2m/3)P , it is clear that the above goes to zero as m ^ oo; uniformly for 
< r < 1. The bound corresponding to (j5.22p with Am replaced with A„ goes to zero at least as 
fast. 

Finally, we sum over Vs. These sets extend over a large fraction of A^, and therefore, they must 
correspond to terms with small norms. Indeed, 



r Y II [^AAZ,r)-^A„SZ,r),a^-iA)] \\dr 

;2s\\A\\ Y E E f^^P II^A„(^,r)||+ sup \\^A^{Z,r)\\) 

.^A ,.,,,^7^ V0<r<l 0<r<l J 



xeA^/ayeA- ZcF 
' x,y£Z 

<4sP||||vI/|| E E F^idi^^y))- 

3^'6A,„/3 ?/eA^^/3 

As is proven in Theorem 14.81 the function F^JJ{r) = u^{r)F{r) for some > and r large enough. 
Thus the sum 



(5.23) 



E E F^{d{x,y)) <\\F\\\Am/3K{bim/3) 
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which goes to zero as m ^ oo. We have shown that all terms vanish in the limit, and therefore, 
the sequence {a^"{A))n is Cauchy as claimed. □ 

The above result establishes the existence of the spectral flow in the thermodynamic limit, and 
we have denoted that limiting flow by . Arguments similar to those used in the proof of Theorem 
15.21 show that is also the thermodynamic limit of the flows generated by the interaction ^'r(s), 
defined in (I5.16p . restricted to the sequence of finite volumes A.^. This is not a surprise since, as 
the next proposition shows, ^r('S) is the limit of ^a(s) as A — ?> T. In this proposition, we consider 
the interactions ^'a('S) as functions defined on the power set of A, ^(A), with values in the algebra 
of observables. As such, we can consider the interactions obtained by restriction to a subset of 
P(A), such as ^'a(s)|p(Ao), for Aq C A. 

Proposition 5.3. For any finite A C T and Z C A, the following estimate holds 



(5.24) \\^AiZ,s) - ^riZ,s)\\ < C\\d<^\\a\Z\^K{diZ,A-)) 

where 



C = 4F(0) J\\W,\\, + Y,VG^) 



n>0 



Let (Am)^ he a sequence of finite volumes satisfying the properties \5.4\) - Then, for any j3 € (0, 1), 
one has 

(5-25) lim ||^'A^jp(A p)-^t\v{K ^)||f*=0 

Proof. To prove the estimate ()5.24p for fixed Z, we apply Lemma l5.ll with A = 4*'^(s) and then 
Lemma l4.1l as follows: 

\\^K{Z,s)-^>r{Z,s)\\ < E ||AX(<i>'y(.),.)-Ap(ci>'^(s),s)|| 

Y,n>0: 
Yr,=Z 

< 4 \Y\WY{s)\WGin - l)K{d{Y,A-)) 




< 4 Vll^7lli + E^^ VK{d{Z,A-)) J2 lYm'yis) 



n>0 / YcZ 



(5.26) < 4||a$||aF(0) ^llW^^lli + E VW) \Z\y'K{d{Z,A-)), 

\ ri>0 / 

which is the claimed result. To prove (|5.25p is now a straightforward application of (j5.24p and the 
properties of the function K defined in Lemma 15. li □ 

Proposition 5.4. The spectral flow for the infinite system has the following properties: 

i. (af)^gjQ is a strongly continuous cocycle of automorphisms of the C* -algebra of quasi-local 
observables, and it is the thermodynamic limit of the finite-volume cocycles generated by the 
interaction ^r(s)- 

ii. oFg satisfies the Lieb-Robinson bound 



(5.27) II [af (^),5] II < 2p||||S||min 



l,g{s) E F'^id{x,y)) 

xeX,yeY 
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for any A € Ax, B G Ay, and < s < 1, with g given by 

^2\mcF^\t\ otherwise. 



(5.28) Cf^ ■ g{t) 



and the quantities Fq,, Cp^, and as given in Theorem \4.5\ 

iii. If P is a local symmetry of^, i.e. , an automorphism such that /3{^{X, s)) = ^{X,s), for 
all X CF and s E [0, 1], then (3 is also a symmetry of a^^, i.e. , ct^o^ = oF^ for all s € [0, 1]. 

iv. Suppose T is a lattice with a group of translations {Tx),j. and {ttTx)x representation of 
the translations as automorphisms of the quasi-local algebra Ar- Then, if ^ is translation 
invariant, i.e. , ^{Tx{X), s) = 7rT^(<I>(X, s)), for all X C T, and s G [0,1], then 
commutes with ttt^, for all x and s. 

Proof. All these properties follow from the preceding results. □ 

5.2. Automorphic equivalence of gapped ground states. We can now describe more precisely 
the problem of equivalence of quantum phases discussed in the introduction. Let Sa{s) denote the 
set of states of the system in volume A that are mixtures of eigenstates with energy in I{s) and let 
iS(s) be the set of weak-* limit points as n — )■ oo of S\^{s). Note that these sets are non-empty. 
The result of Section [2] immediately implies 

(5.29) 5A„(s)=5A„(0)oa^, 

where a^" is the automorphism defined in (14. lip . In Section [4] we proved that a^" satisfy a Lieb- 
Robinson bound with a uniformly bounded Lieb- Robinson velocity and decay rate outside the 'light 
cone'. In the previous subsection we obtained the thermodynamic limit of these automorphisms 
leading to the cocycle which automatically satisfies a Lieb-Robinson bound with the same 
estimates for the velocity and the decay. The following theorem states that (I5.29P carries over to 
the thermodynamic limit. 

Theorem 5.5. The states uj{s) G S{s) in the thermodynamic limit are automorphically equivalent 
to the states a;(0) G 5(0) for all s G [0, 1]. Indeed, 

(5.30) S{s) = 5(0) o of 

Moreover, the connecting automorphisms can he generated by a s- dependent quasi-local inter- 
action ^{s) with \\^\\p^ < oo, where the norm is defined in ^5. 5\) . then satisfies the same 
Lieb-Robinson bound as in Theorem \4.5l 

Proof. This is a direct consequence of (j5.29p . theorem 15.21 and the lemma below. □ 

Lemma 5.6. Let {crn)n be a strongly convergent sequence of automorphisms of a C* -algebra A, 
converging to a and let {oJn)n be a sequence of states on A. Then the following are equivalent: 

i. Un converges to uj in the weak-* topology; 

ii. iOn° cr converges to uj o a in the weak-* topology; 

iii. ujn ° converges to uj o a in the weak-* topology. 

Proof. follows immediately from the fact that a and are automorphisms. Now if (ii) 

holds, the second term of 

|(w„ o an){A) - (wo a){A)\ < ManiA) - a(A))| + \uJn{(T{A)) - u}{a{A))\ , 

vanishes. So does the first one 

\uJn{(Jn{A) - a{A))\ < WunWhniA) - a{A)\\ 

since ujn are states, and therefore (iii) holds. A similar argument yields (iii)=>(ii). □ 
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In the recent literature [IOl[TT], a 'ground state phase' has been defined as an equivalence class 
of ground states with the equivalence defined as follows: the states loq and uji are equivalent {i.e. , 
belong to the same phase) if there exists a continuous family of Hamiltonians H{s), < s < 1, such 
that for each s, H{s) has a gap above the ground state and ujq and uji are ground states of H{0) 
and H{1), respectively. As an alternative definition the authors of [10] state that ujq and ui should 
be related by a 'local unitary transformation'. With Theorem 15.51 we provide precise conditions 
under which the first property implies the second. At the same time we have clarified the role of 
the thermodynamic limit left implicit in the cited works. 

Based on Theorem 15.51 it seems reasonable to define the ground states of two interactions <&(0) 
and ^{1) to be in the same phase if there exists a differentiable interpolating family of interactions 
< s < 1, such that there exists a > for which \\^\\a + ||c^*J'||a < co, and if the spectral gap 
above the ground states of the corresponding finite- volume Hamiltonians H\^{s) have a uniform 
lower bound 7 > 0. The increasing sequence of finite volumes Am should satisfy a condition of 
the type ()5.4I) . One should allow for a space of nearly degenerate eigenstates of which, in 
the thermodynamic limit, converge to a set of ground states S{s). We have proved that under 
these conditions the sets of thermodynamic limits of ground states are connected by a fiow of 
automorphisms generated by a quasi-local interaction with almost exponential decay and satisfying 
a Lieb-Robinson bound. We believe that these are sufficient conditions for belonging to the same 
gapped ground state phase. More work is needed to identify necessary conditions. 

We remark that a 'ground state phase' should be defined as an equivalence relation on simplices 
of states of a quantum lattice system. This is an equivalence of sets of states rather than of models 
because it is possible that different quantum phases coexist as ground states of one model, while 
the same states also appear as unique ground states of other models. Examples of this situation 
can easily be constructed using frustration free models in one dimension with finitely correlated 
ground states, also known as matrix product states [171 ES]- In particular, if 5(s) denotes the 
set of infinite-volume ground states of a model with parameter s, the relation S{s) = 5(0) o a^, 
does not imply that the states in the sets <S(s) are automorphically equivalent among themselves. 
E.g., if for a model with a discrete symmetry we find that symmetry broken states coexists with 
symmetric states, cannot map these two classes into each other. In general, as emphasized in 
Proposition 15.41 the Og we constructed posses all symmetries of the Hamiltonians. 

There are plenty of examples of models to which our results apply. Clearly, the various pertur- 
bation results mentioned in the introduction provide many interesting examples of sets of models 
with ground states in a variety of types of gapped phases. Another class of examples is provided 
by the rich class of gapped quantum spin chains with matrix product ground states. In Yarotsky's 
work [51j it is shown how perturbation theory around a matrix product ground states can be ap- 
plied to connect these two classes of examples. Exactly solvable models with gapped ground states 
depending on a parameter, such as the anisotropic XY chain [33], is another set of examples. More 
recently, stability under small perturbations of the interaction was proved for a class of models 
with topologically ordered ground states [7]; these include e.g. Kitaev's toric code model [30]. Our 
results are also applicable to this class of models. It seems likely that other applications will be 
found. As an example of an application left to be explored, we mention that the existence of a 
connecting automorphism of the type can provide a means to distinguish true quantum phase 
transitions from isolated critical (i.e., gapless) points around which it is possible to circumnavigate 
with suitably chosen perturbations. 
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